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16. ABSTRACT 

During the first quarter, progress was made in all ttirce areas which were included 
in the scope of work. The following report describes the work done in three sections. 

;\ sputtering system was developed to deposit aluminum and aluminum alloys by the 
D.C. sputtering technique. This system is designed for a high level of cleanlinc.ss 
and for monitoring the deposition parameters during film preparation. This system is 
now ready for studying tiie deposition and annealing p.irameters upon double- level metal 
preparation. A new teehnique recently applied for sem iconducti'r analysis, the finite- 
element method, was studied for use in the computer modeling of two-dimensional MOS 
transistor structures, and it was concluded lliat tlie mettuul lias not been suff icient Iv 
well developed for confident use .u thi.s t imt>. An algoritiim lia.s liei-n tleve It'ped fei' 
implementing a computer study wliieh is based upon the f i ni t e-d i 1 f erenec method. The 
program which was develtiped during the la.st scope of wi.>rk was muiiil ied and used ti' 
calculate redistribution data for boron and pltosphorous wliich iiad bci'n predoposiied 
i\v ion-implantation with range and -straggle conditions typical of tliose used a M.RFG. 
Data was generated for -111' oriented films witli red 1st r i hut i<m in N.^, dry 0^ .and 

sto.am ambients. D.ata will Ih* giviMi in the I inal repi'rt showing both t w(' t1 inu'ns ion.i ! 
effects and the evolution ot the junct ion depth, .sheet resistance and integr.iteil do.se 
with redistribution time. 
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1. DOUBLE-LEVEL METALLIZATION TECHNIQUES 


Moat of the work during the first quarter was directed toward 
the development of a sputtering system for preparing aluminum and 
aluminum-alloy films. A photograph of the completed system is 
shown in Figure 1.1. Briefly, the system consists of the following. 

The basic vacuum system is a Varian-NRC 6 inch oil diffusion 
and mechanical roughing pump equipped for automatic and manual 
operation. The diffusion pump is equipped for LN^ cooling of the 
cold trap. The sputter gun and power supply were obtained from 
Sloan Technology. The sputtering chamber was designed and built 
at Mississippi State. It consists of a Corning 12x18 glass cylinder 
and an aluminum top plate machined to accomodate up to three sputter 
guns. A cold cathode discharge gauge was constructed and installed 
in the baseplate of the chamber to measure the pressure during the 
sputtering operation. A throttle valve with several threaded holes 
for accomodating plugs is operated by one of two mechanical feed- 
throughs in the base plate. A lift mechanism with a reversible 
motor was designed and contructed for raising the top plate and 
sputter guns. 

Inside the chamber is equipped with a rotating table which 
accomodates up to eight wafers of l*s-2" in diameter. The table is 
driven by a vacuum sealed shaded pole motor through a magnetic 
coupling at 7 rpm. The entire motor-tabl<“ assembly is rotated by 
a chaln-sprocket-mechanlcal feedthrough arrangement through three 
sputter-gun and two mask positions. The mask is attached to the 
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rotating assembly and provides one hole through i»hlch the sputter- 
gun deposits metal on the wafers. A crystal film thickness sensor 
Is located beneath the spiittcr-gun and receives a deposit through 
a hole at an unused wafer position on the table. 

The sputtering chamber opens into a class 100 clean bench in 
order to maintain a high level of cleanliness. The system is 
located in the metallizing and bonding room of the microelectronics 
laboratories in Simrall Engineering Building, and this room was 
designed with air conditioning and filtering units to maintain a 
class 10,000 environment. The entire system has been constructed 
and checked out and is ready for depositing sputtered films. 

The installation of a six-tube Thermco Ranger diffusion furnace 
was completed with the addition of a venting system for exhausting 
the scavenger boxes. All that remains to be done is to line the 
tube and connect the nitrogen ambient source in order to anneal 
the aluminum films. 

2. TWO-DIMENSIONAL MODELS FOR MOS TRANSISTORS 

The work done during the first quarter was directed toward 
preliminary Investigations of numerical schemes and computational 
algorithms for solving the semiconduction equations for a two- 
dimensional field. 

A recent report lins described the application of the flnlte- 
elenicMJt metliod Lo tlie analysis of a JFET * . Tin- f ini te-fleinent 


method has been used for some Lime in salving problems in mechanics 
and elasticity; however, it has only recently been applied to 
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••mlconductlon problems. This method has the power to treat some 
problms, such as eigen-value problems, for which the finite- 
difference nwthod is awkward if at all applicable. It can also 
be applied to the solution of field distributions governed by 
partial differential equations, and one of the most attractive 
features as compared to the finite-element method is purported to 
be the ease of treating non- rectangular geometries and Irregular 
boundaries. For example, the geometry of the VMOS structure could 
be accomodated. It was decided to further Investigate this 
technique. 

In order to better understand the applicability of the method, 
it was applied to a one-dimensional linear diffusion problem. This 
simple problem is one for which familiar results are available for 
comparison and at the same time taxes the finite-element method. 

Ii> its most valid form, the finite element method is applicable to 
variational problems in which a true minimum of an energy-related 
function exists. Such a minimum does not apply for the semi- 
conductor problem in which current flow occurs by diffusive and 

conductive mechanisms. It has been proposed that a "weak form", 

2 

the so-called Galerkin method, be applied to such problems. The 
typcial seminconductor problem is a non-linear boundary value 
and initial condition problem of which the linear diffusion 
problem Is a very special case. In the example chosen, the dif- 
fusion variable, u, obeys; 
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u(x»t)lj^Q - u^ - 1 

(•) 

(2.1) 

“*<*•'> lx.« ■ « 

(b) 


■ 1 » x-0 
0 , x>0 

(c) 


u^ “ u ■ f (x,t) ■ 0 

t XX 

0<x<a 

(2.2) 


The Gelerkin foraulatlon of chit problra ie: 

Jq ( u^v - - fv ) dx - 0, (2.3) 

%dtere v(x.t) represent a *'crlal function" which is used to 
approximate u(x,t). The finite element methods uses a set of 
"hill functioiis" as illustrated in Figure 2.1 to construct the 
v(x»t) approximation. Two of the popular hill functions are 
the Hermite bicubic and the bilinear functions which are illus- 
trated in the figure and were used in the example. The final 
form of the approximate solution is; 

N 

v(x,t) ■ E qi(t) ^,(x) . (2.4) 

i-1 ‘ 

On the node points the solution is approximated by the set {qi(t)} 
for the type of hill functions which overlap as illustrated in 
Figure 2.1. Solution for the set (q£(t)} is then analogous to 
solving for the set {uj^(t)} on the node points using the finite 
difference technique. The equations for the set (q|(t)} are 
obtained by substituting (2.4) into (2.3). 
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(a) Btilnaar hill functions 



« - ■ , j ■ k X 

0 12 3 4 

<b) Bicubic 'Hsfiilta) hill functions 


Figure 2.1. Illustration of hill functions used 
In finita aiaoMt oathod. 
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.f j) 


dx 


.1 - 1.2,3,. .....N (2.5) 

From (2.5) O act of tine differential equations la obtained 
which Is solved uaing an implicit ntmerical method. 

The solution of the problem posed by the example la closely 
approximated by the erfc function in the range 0 £ x < 3 if 
a ■ 6, and this solution was used to compare the accuracy of the 
finite element and finite difference methods. Figures 2. 2-2.4 
show the max laws error as a function of the reciprocal of the 
number of grid points and the size of the tine step, i.e. 

■ dt/Ax^ . 

The error obtained in the solution by the bilinear finite 
element method is very nearly the same as that obtained with the 
finite difference method. This was not surprising because the 
system of equations for iq^(t)) and (u^(t)) were quite similar. 

What was surprising was that the Hermite bicubic finite element 
produced such poor results, although this surprise was based upon 
the intuition that since this element was more difficult to use 
it should provide some reward for the difficulty. 

This experience resulted in some skepticism that the finite 
element siethod would be effectl'^-' for semiconductor problems. A 
XBore recent paper has reinforced this attitude.^ This paper Indi- 


cates that the proper formilations of the semiconductor problem 
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Figure 2.2. Error for finitoediffortnct vs, roclproctl of toUl 
moiber of grid points. 
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Error for finite-bicubic-element ys, reciprocal of 
total number of grid points. 


1 / TOTAL GRID POUTS 


for £h« fiaito eioMat approach raBsin to ha donnitratad, and, 
in agreanont with our oboervatlon«» point <nit that tha application 
of 6alarkin*t nathod la aubjact to akaptlciaai, Tharafora, it was 
concluded that the further work should be based upon the finite- 
difference method which wa have used before although the finite- 
element method is intriguing and may be further developed in the 
future. 

The second phase of this program is underway to develop a 
computational program with which to generate data from a model. 
During the past quarter* the major emphasis was upon deriving 
a simple two-dimensional algorithm which could be econcsBlcally 
used In a simulation program during its development. It is 
plausible that this algorithm can be later modified to become 
adequate for modeling of other effects which are significant In 
certain situations, e.g., short channel effects, avalance break- 
down, etc. 

During the past quarter an algorithm has been developed based 
upon the usual assumption chat the mobile carriers are Included In 
an Infinitesimally thick layer of charge at the S1-S102 interface. 
The current flow Is then described by a one-dlmenslonal equation; 

where 1 and S are channel current and charge per unit channel width, 
E Is the tangential Interface field and and p. are the mobility 

8 M is 

factors accounting respectively for gate modulation and hot 
electron effects. This equation Is solved Iteratively with 
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Poisson's equation which includes two dinenslonal effects: 

Vhm -p/e . (2.7) 

Equation (2.7) will be solved for one segment of a periodic struc- 
ture with respect to x and equation (2.6) will be integrated to 
produce auxiliary equations for the boundary conditions at the 
SI-SIO 2 Interface. It will be assumed initially that S ■ 0 at 
X “ a, the point at which the normal component of the Interface 
field changes sign. The solution algorithm then proceeds in an 
iterative fashion to solve (2.6) and (2.7) simultaneously for the 
potential distribution from which the current, i, is untlmately 
calculated as a function of the gate, drain and body voltage. 

The model at this point admittedly has some short-coming, 
mainly due to the neglectlon of generation-recombination mechanisms. 
Therefore, it will not treat impact avalance and bulk generated 
leakage currents. It is believed at this time that such effects 
can be treated by adding another iterative loop to the algorithm. 

The major emphasis at this point will be the development of a 
program for input and output data management and including sub- 
routines which generate data Internally within the program and 
solve svstems of equations which will be encountered in imple- 
menting the algorithm. 
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3. REOISTRIBUTKMI DIFFUSIONS F(» ION- IMPLANTED PREDEPOSITS OF 
BORON AND PHOSPHORUS IN SOS FILMS. 

The objective of this work was to produce curves describing 
the variation with diffusion time and temperature of the junction 
depth, sheet resistance and integrated Impurity dose. This data 
has been generated for boron and phosphorus redistributed in 
nitrogen, dry oxygen and steam ambients for <111> oriented SOS 
films. The following section presents discussions of the implan- 
tation and redistribution model, further program develop, the 
computational procedure and of the computed results. 


3.1 The Redii rlbution Model: 


There are three aspects of the redistribution model which are 
considered: (a) the Implanted profile, (b) the oxidation model, 

and (c) the dlffusivity model. 

(a) The Implanted profile. The implanted profile is 
described by the gausslan function, 

V R 2 

c<y) = f'w Pn ! . <3.i) 

max ( A Kp ' 

where C is the concentration, y is the distance from the entrant 

silicon surface, R is the range and AR is the straggle for the 

P P 

implant. The peak concentration, is related to the implant 

'’imp 


c 

max 



// 2ti ' AR 

P 


(3.2) 


Redistribution data has been generated for the following ondi t Ions 
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Q* j 5x10^^ , 10^^ , 5x10^^ , 10^^ cm“^ 

imp 

R : 0.2735 Mm 80 keV boron implant. 

P 

0.1727 " 150 keV phosphorus. 

AR : 0.0665 pm 80 keV boron Implant. 

P 

0.0440 " 150 keV phosphorus. 

The doses are light to nnderate resulting in concentrations no 
18 —3 

heavier than 6x10 cm , and the range-straggle values are typical 
of those employed at MSFC. It is assumed that all of the ions 
become activated shortly after redistribution begins and thereby 
diffuse by a substitutional mechanism involving vacancies. 

(b) Oxidation Model: The oxidation model is assumed to be 

the same as for bulk silicon and the data of Deal et.al. has been 
used to calculate the oxidation rate according to: 

dx 

-7^ - B/ (2x + B/C) . (3.3) 

dt o 

where B and C follow Boltzmann- like temperature dependences. 

Figures (3.1) and (3.2) illustrate the oxide thickness dependence 
upon time and temperature for both dry 0^ and steam ambients. 

During the oxidation, the silicon film thickness is reduced 
according to: 

W = W - ax , (3.4) 

o o 

where W is the initial film thickness, taken to be 1 pm, and 
o 

a ■ 0,45 is the ratio of the densities of Si 02 to silicon. 
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Redistribution data is given for W ■ 1 mid and an initial oxide 

e 

thickness of x ■ 300 A. 
o 

(c) Diffusivity Model: The diffusivlty model for boron was 

discussed in an earlier report^ and it Includes a linear depen- 
dence of the diffusivlty upon the vacancy concentration as well 
as the field-enhancement effect. The diffusivity model for phos- 
phorus includes only the field-enhancement effect which is suf- 
ficient to describe the non-linear behavior of phosphorus diffusions 

19 -3 7 

at concentrations lower than 10 cm as shown by Barry and 

g 

Fair and Tsai . The dlffuslvlty-temperature dependence is after 
Falr^ and Fair and Tsai^ adaptation of data by Ghostagore^^, For 
either boron or phosphorous the effective diffusivity is given by: 


where, 


and, 


Dgff “ D(u) X (1 + u//„2 + 1 ) , 

u » C/ 2 n^, 

D(u) = 0* u , for boron, 

D 

= D* , for phosphorus. 


(3.4) 


(3.5) 


(3.6) 


and where n^ is the intrinsic carrier concentration at the diffu- 
sion temperature and D| and D* are the intrinsic dlffusivities of 
boron and phosphorus: 


“b 


3.17 exp ( -3.59eV/ k^T) cm /sec., 


D* 

P 


3.85 exp (-3.66eV/ k_T) 

o 


(3.7) 
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3.2 Further Program Development 

The program which was used to generate the data has been 
described In detail in an earlier report. It was noted that the 
program was developed In such a way that one could take advantage 
of a normalisation procedure for predeposition diffusions and 
generate data applicable to different film thicknesses. However, 

It is not possible to gain such an advantage for redistribution 
diffusions involving ion- Implants or growth of an oxide. Then 
the program was used to generate data. It was discovered that some 
other features of the program are extraneous unless further 
refined. 

The program was developed to account for both thin and thick 
oxides such as would be encountered in some practical situations. 
However, such a simulation requires the Incorporation of a 
warped grid system, a modification which would require considerably 
more effort. Therefore, the variable oxide feature is of limited 
value at this time, since the program, at best, only approximates 
Che conditions for growth of a very thin oxide during redistribution. 

A modification was made which allows accurate treatment of 
redistribution under oxidizing conditions when only a single oxide 
thickness is involed. The original program treated the oxidation 
process with regard to Che boundary conditions; however, unlike 
the case of bulk silicon, one must also account for the reduction 
of Che silicon film thickness. This feature is now included in 
the program. During the simulation of a redistribution in an 
oxidizing ambient, the vertical grid spacing continuously shrinks 
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while the horizontal grid spacing la constant. The amiificatlon 
does not show up on logic flow diagrams at the level of detail 
which has previously been given. For ccanpleteness, a new listing 
of the affected main and sub-programs will be given In the 
Final Report. 

3.3 Computational Procedure; 

The program described in an earlier report, and modified 
as described in the preceding, was used to generate the data. Ibo- 
dlmenslonal data was obtained in the form of Isoconcentration 
contours for typical situations. The bulk of the data which can 
be correlated with experimental measurements is generated using a 
qtiasi-one-dlroensional model in a manner described in a previous 

g 

report . A brief review of the procedure is given in the following. 

For generation of sheet resistance. Junction depth and 
integrated impurity dose data as a function of time and temperature, 
only a one-dimensional profile need be calculated. This is 
accomplished by making the horizontal grid only three units wide 
but keeping the field six film thicknesses wide. Periodic 
boundary conditions for the horizontal dimensions are employed in 
the program and result In a calculation which produces the vertical 
profile equivalent to a none-dimensional model. Thus without 
sacrificing the generality of the program for treating two- 
dimensional cases, the amount of computing time is drastically 
reduced when the data that is desired does not require the full 
power of the program. 
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The vorticnl grid vailos ffom thirty one to sixty one 
points as required for accuracy In details of the proiile, and 
most of the data is not sensitive to the number of grid points 
used if the number is chosen in this range. For the purpose of 
illustrating the unusual nature of phosphorous profiles, the 
larger number of points was required. 

3.4 Discussion of Results 

First, some of the unusual behavior of redistribution dif- 
fusions in SOS films will be discussed in this section. Next, the 
format for the calculated curves will be discussed. 

Figures (3.1) and (3.2) illustrate the oxide thickness growth 
and silicon film thickness reduction as functions of time for 
<111> silicon slims oxidized in steam and dry O 2 ambients. The 

ft 

evolution traced beginning with an initial oxide of 300 A thickness 
on an SOS film of I »im initial thickness. The curves arc shown 
for four temperatures. The data ore necessary for Interpreting 
some of the results for simulated redistributions. 

Figures (3.3) and (3.4) show impurity profiles for boron 
and phosphorus Implants being redistributed in a steam ambient at 
1000 deg. C. The profiles are all plotted with a coomon origin 
as would be the case for experimentally derived profiles where 
Sl-SlO^ interface would serve as the logical origin. However, 
the profiles are normalized with respect to the film thickness 
which 1s of course shrinking. The boron profiles are not unusual 
but show the weliknown leaching effect due to segregation into 
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the oxide. The phosphorus profiles show the effect of impurities 
being rejected from the oxide. There is a pile-up of impurities 
in front of the advancing S1-S102 Interface and then a dip which 
eventually disappears. It is easy for one to draw an erroneous 
conclusion from observing the profiles, because It appears that the 
integrated dose should Increase or at least remain constant and 
the sheet resistance should decrease with time. This Is not true. 
Although the segregation coefficient favors phosphorus in silicon 
vs. S 102 , eventually all of the phosphorus will be in the Si 02 
when the SOS film is corq>letely oxidized since the model assumes 
that there is no diffusion into the sapphire. 

Figures (3. 5-3. 7) Illustrate the behavior of the junction 

migration, sheet resistance variation, and integrated impurity 

dose variation over a long period of time. All of the curves are 

plotted with respect to normalized time, and true time is obtained 

by multiplvlng by the normalizing time value given on the plot. 

Junction depths are in microns, sheet resistance values are in 

_2 

ohms, and dose values are in cm units unless otherwise marked. 

The curves are given in the typical format for all of the data. 

For an ion-implanted profile, there are in fact two junctions 
until one of the junctions emerges at the Sl-Si02 interface. There- 
fore, the sheet resistance values are for the buried layer until 
the front junctlo.. disappears. This typically occurs in a short 
time compared with that for through-diffusion of the back junction. 
Figure (3.5) illustrates the through-diffusion of the back junction 
for the heavier doses. This always .occurs for redistribution in 
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a nitrogen ambient but not necessarily so for an oxidizing 
^bient. After the through diffusion, or even before for light 
doses, the junction depth will eventually decrease due to the 
reduction of the film thickness or due to the relatively slow 
advancement of the junction with respect to the moving S1-S102 
interface. 


ytV of 
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